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SUMMARY 


A closed-form asymptotic expression for the mutual admittance between 
two circular waveguide-fed apertures is developed for large element 
spacings. A comparison with calculations obtained by a numerical 
integration method indicates that the present analysis gives good results 
for spacings greater than two wavelengths . 

INTRODUCTION 

In the analysis of mutual coupling effects between radiating elements 

of a finite size antenna array, a calculation of the complex mutual 

admittance (or impedance) between each pair combination of the array is 

2 

needed. For example, an array of N elements requires N values of 
complex mutual admittance to form the complex admittance matrix from 
which the complex scattering matrix is obtained (see ref. 1). 

A considerable reduction in the actual number of calculations can 
be achieved by taking advantage of the periodicity and symmetries of the 
array configuration; however, there still exists the need to calculate 
the mutual admittance between elements which are many wavelengths apart. 

An earlier analysis (ref. 1) requires the numerical evaluation of an 
integral for the calculation of mutual admittance between elements of the 
array. It has been found that the computer time required for the numerical 
integration increases as the distance between the pair of elements 
increases. The purpose of the present paper is to present an asymptotic 
evaluation of the integral for large element separations. 
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SYMBOLS 

radius of aperture 

\ 

electric field vector 

E^U,y) 

functional form of 


Fourier transform of E^{x,y) 

\ >\ 

X and y components of E^ 

X y 

E (a,3),E (a,3) 

Fourier transforms of E^ end E^ 

X y 

■^x V 

f(3) 

functional notation 

f'(3) 

first derivative of f(3) vith respect to 3 

f"(3) 

second derivative of f(3) with respect to 3 

G^(B) 

quantity defined "by equation (T) 

GgCP) 

quantity defined by equation (8) 

g(3) 

functional notation 

^2 

magnetic field vector 

EgCx.y) 

functional form of 

H„(k ,k ) 
2 X* y 

Fourier transform of ^^{xjy) 

Hg ,H 

X and y components of Hg 

X y 

(a,3),H2 (a,3) 

Fourier transforms of Hg and Hg 

X y 

X y 


Bessel function of the first kind of order n ; 

J^(v) 

first derivative of with respect to w 


k = 2 tt/X 
o 

plane wave pi'opagation constant in free space 

k jk ,k 
x’ y’ z 

wave propagation constants in x,y,z direction 


vnit outward normal to aperture 1 

R 

center-to-center spacing of apertures 


f 


i 






r,0,4> 


S (k 3R) 

U O 

S2(k^3R) 

x,y,z 



z 

a 


3 



spherical coordinate variables 
area of aperture 1 
q,uentity defined by equation (ll) 
quantity defined by equation (12) 
modal voltages 

variables in Cartesian coordinate system 
complex mutual admittance 
unit vector in z direction 

angular Fourier transform variable in cylindrical 
coordinate system 

normalized radial Fourier transform variable in 
cylindrical coordinate system 
infinitesimal value 
permittivity of free space 
free space wave length 
permeability of free space 


THEORY 


The mutual admittance between two waveguide-fed apertures can he 
determined hy the electromagnetic reaction between the equivalent electric 
and magnetic currents (ref. l), i.e. 


’^12 = vX • Si (1) 

1 2 Si 

where Y^^ mutual admittance between apertures 1 and 2, and Vg 

are the normalized modal voltages, is the area of aperture 1, n^ is 
the outward unit normal to S^, E^ is the assumed electric field excitation 
of apertiire 1, and Hg is the short circuit magnetic field due to an assumed 
excitation of apert\ire 2. 

If the two apertures lie in the x,y plane, ets shown in figures 1 and 2, 
then the magnetic field outside aperture 1 due to the excitation of 
aperture 2 can be expressed in the form 

“ilk z 

Hg = Hg(x,y) e (2) 

where k is the wave propagation constant in the z direction. Then (l) 
can be written as 

• z dx dy (3) 
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V V 
'' 1^2 


// 

S, 


E, 




^'“^Hg(x.y)e 




The evaluation of (3) can be facilitated by letting 6 = r cos 6 and using 
a form of Parseval’s theorem (ref. 2, eq. (C-15)). Then 






\2 = 


lim. 




// 


-jk r cos 0 

Ei(X^,ky) X H2(-k^.-ky) e * 


• z dk dk 
X y 




where k^ and k^ are the wave propagation constants in the x and y 
directions. 

The mutual admittance can also he expressed in terms of the two- 
dimensional Fourier transforms of the assumed aperture electric fields 
and the solutions to the wave equations for the region outside the aperture 
plane. The details are outlined in ref. 1; therefore (U) becomes 


k I— 


\2 = 


lim 


“ 2 tt -jk 

/ / 


..v^T 


3*^ cos 6 


3=0a=0 

- 7 ■ /Et (a,3)cos a + E, (a,3)sin a)(E„ (a,-3)cos a + E„ (a,-p)sin a) 

V /V^x / 

+ Jl - ' 0 ^ |E^ (a,3)sin a - E^ (a,3)cos ^(^2 (^>”^)sin a - (a,-6)cos aj 

\ X y / \ X "^y / 


• & d0 da 


(5) 


where a change of variables has been made in the transform domain such that 

k = k 0 cos a and k = k g sin a. 

X o y o 

If both apertures are excited by the circular waveguide TE^^ mode, 
then the mutual admittance becomes (details are in ref. l) 


Y =2 
12 ^ \l\l 


lim f 
8-^/2 8=0 


-jk r Hi - 3^ cos 6 


L(l.8Ul)" - Ij 

Gj^(B)(jQ(k^3R) + J2(k^3R)cos 2<j>] 

+ G^(3)(j^(k^3R) - Jg{k^3R)cos 2<f)) 


3 d3 


( 6 ) 


where 


\(3) = 1=^=^ 
^ J±- 3^ 


J^(k^a3) 


(T) 


,(3) = 


j;(k a3) 

k a 1 o 
o 

a.84i\"2 p2 


L(W- 


( 8 ) 


Now Tising the semi convergent series of Hankel (see pages 137 end 138 of 
ref. 3), Jq( 1!:^3R) and J^(k^3B) can be expressed as 


2' o'' 
r 


J (k 3R) = e 
° ° ^ 


1 IT 


jk^3R 

S S^(k SR) + j 

vor ° ° 


-jk^SR 

S (-k SR) 

/Or ° ° 


-.1 ^ 


J^d^^BE) = — e 

/27T 


jk^SR 

^ S„(k BE) + ,1 S„(-k, BE) 

,Or ^ ° /5Tbe ^ ° 

o o 


(9) 


(10) 


where 


S^(k^BE) 


Satk^BH) 


1 ^ \ 

_ -2_ 1 

( ^ V 


128 ’ 

Uo®"'/ 


Jl. 

1024 


(V") * ■ 

^ 8 ^k^SRj 128 \^k^3R/ ^ 1024 i^k^SRj 


(11) 

( 12 ) 
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The mutual admittance then hecomes 


1 m 

jr~ ^ 4 W rjl - cos 0 

p j_2. 3, e lim f e ° ’ 

'1^ (l.84l)^- 1 o 


'o 1 


Jlc^SR 
e 

I 


( g ^{&) - G^(3))s2(k^0R)cos 2<j. + (GgO) + G^(3))s^(k^0H) 




(g 2(3) - G,(e))s2(~k^6R)cos 2(f> + + G^(3))s^(-k^0R) | 


-jk^3R 

If 3 is substituted for -3 i^i terms which contain e and ve let 


R ~ T sin 0, (13) can he expressed as 


Y .= p L° _i_ _® r e 

(l.84l)^- 1 


-jr^k^3 sin 0 - k^/l - 3^ cos sj 


* ~ 

vk 3r sin 0 ^ 

o 

+ (Gg(3) + G^(3))s^(k^3r sin 0) 3 d3 ( 

An integral of the form in (l4 ) can he evaluated asymptotically for 
large values of r (ref. 4), i.e. 


~ Jrf(3) jrf(S^) 

/ g(3)e d3 ~ e g(3 ) 

-00 *' 0 


f'(e) 


= 0 


3=3, 


(16) 


f"(3^) = f"(3) 


i- 0 


3=3, 


From (15) > we have 


f{3) = ^ ~ \/l - 3^ cos 0 j 


(IT) 


(18) 


therefore. 


3 5= -spin 0 

o 


(19) 


f(3 ) = - k 


( 20 ) 


= -t; 

cos 0 


( 21 ) 


and the asymptotic evaluation of the mutual admittance for large spacing hecomes 
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J-(h a) 
1 o 

2 

“ -jh R 
e ° 

IVo 

(1.81+1)^ - 1 



^o» 


2 

2Bin24,-^(3-35=os24,)(jij) 

\ o / _ 

f ' ' r 

Hv) i * 2*) - tMiT (5 - 21 2'1')(st) 


( 22 ) 
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RESULTS 


Calculations are given in figures 3 and 4 for the mutual admittance 
of circular waveguide-fed apertures in both E-plane (4> = ir/2) and H-plane 
(<|i = O) coupling configiirations. A comparison is made between the 
asymptotic evaluation (equation 22) and the n\mierical evaluation (ref. 1) 
as a function of aperture spacing. Good agreement with the mutual 
admittance obtained by numerical integration is indicated for element 
spacings greater than two wavelengths. 

CONCLUSION 

The present analysis has been shown to give good results for elements 
spaced greater than two wavelengths. The judicious use of the asymptotic 
expression for mutual admittance together with the numerical integration 
method covld result in a considerable savings in computer time •vdaen 
calculating the effects of mutual coupling in a large planar array. 


( • 
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Figure 1. Circular waveguide- fed apertures in x,y plane. 



Figure 2. Cross sectional drawing of circular waveguide-fed apertures. 




Figure 5. Calculated mutual admittance of two circular waveguide-fed 
apertures for E-plane coupling. 
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